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Abstract 

. An algorithm for solving first order ODEs, by systematically determining symmetries of the form 

C<") ' [£ = F(x), n — P(x) y + Q(x)] - where £ d/dx + n d/dy is the symmetry generator - is presented. To these 

linear symmetries one can associate an ODE class which embraces all first order ODEs mappable into 
separable through linear transformations {t = f(x), u = p(x) y + q(x)}. This single ODE class includes 
as members, for instance, 78% of the 552 solvable first order examples of Kamke's book. Concerning 
the solving of this class, a restriction on the algorithm being presented exists only in the case of Riccati 
type ODEs, for which linear symmetries always exist but the algorithm will succeed in finding them only 
partially. 



1 Introduction 
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One of the most attractive aspects of Lie's method of symmetries is its generality: roughly speaking, all 
solving methods for DEs can be correlated to particular forms of the symmetry generators pi However, 
for first order ODEs, Lie's method seems to be, in principle, not as useful as in the higher order case. The 
problem is that the determining PDE - whose solution gives the infinitesimals of the symmetry group - has 
the original first order ODE in its characteristic strip. Hence, finding these infinitesimals requires solving 
the original ODE, which in turn is what we want to solve using these infinitesimals, thus invalidating the 
approach. 

For higher order ODEs, the strategy consists of restricting the cases handled to the universe of ODEs 
having point symmetries, so that the infinitesimals depend on just two variables, and then the determining 
PDE is overdetermined. Although few second or higher order ODEs have point symmetries, and the solving 
of such a PDE system for the infinitesimals may be a major problem in itself the hope is that one will 
be able to solve the system by taking advantage of the fact that it is overdetermined. 

One basic motivation in this approach is also that the finite transformations associated to point symme- 
tries are pointlike and these transformations form a group by themselves - not just with respect to the Lie 
group parameter. Actually, the composition of any two point transformations is also a point transformation. 
Consequently any two point symmetries can be obtained from each other through a point transformation. 
Lie point symmetries can then be used to tackle the ODE class all of whose members can be obtained from 
(are equivalent to) each other through point transformations, and which includes a member we know how 
to solve (missing the dependent variable). 

Such a powerful approach however is not useful in the case of first order ODEs - the subject of this 
paper - for which "point symmetries" are already the most general ones. The alternatives left then, roughly 



speaking, consist of: looking for particular solutions to the determining PDE ||, or restricting the form of 
the infinitesimals trying to emulate what is done in the higher order case, such that the problem can be 
formulated in terms of an overdetermined linear PDE 0]. The question in this latter approach, however, 
is what would be an "appropriate restriction" on the symmetries such that: 

• The related invariant ODE family includes a reasonable variety of non-trivial cases typically arising in 
mathematical physics; 

• The determination of these symmetries, when they exist, can be performed systematically, preferably 
without solving any differential equations; 

• The related finite transformations form a group by themselves - not just with respect to the Lie group 
parameter - so that the method applies to a whole ODE class. 

Bearing this in mind, this paper is concerned with first order ODEs and linear symmetries of the form 

£ = F(x), r, = P(x)y + Q(x) (1) 

where {£, 77} are the infinitesimals, the symmetry generator is £d/dx + rjd/dy and x and y = y(x) are 
respectively the independent and dependent variables. Concerning the arbitrary functions {F, P, Q}, the 
requirements are those implied by the fact that (Q) generates a Lie group of transformations. The linear 
symmetries (Q) have interesting features; for instance, the related finite transformations are also linear, of 
the form 

t = f(x), u=p(x)y + q{x) (2) 

where t and u = u(t) are respectively the new independent and dependent variables, and {/, p, q} are 
arbitrary functions of x. Linear transformations (0) form a group by themselves too - not just with respect 
to the Lie parameter. So, as it happens with point symmetries in the higher order case, any two linear 
symmetries (Q) can be transformed into each other by means of a linear transformation (Q), and hence to 
the symmetries (Q) we can associate an ODE class. Since separable ODEs have symmetries of this form, the 
class of ODEs admitting linear symmetries (Q) actually includes all first order ODEs which can be mapped 
into separable by means of linear transformations (|^) . 

We note that in the particular case of polynomial ODEs, e.g. of Abel type, 

y' = h{x) y 3 + f 2 (x) y 2 + h(x) y + f (x) (3) 

(||) actually defines their respective classes]]. Since a separable Abel ODE can be obtained by just taking 
the coefficients all equal, this means that there are complete Abel ODE classes all of whose members 
can be transformed into separable by means of (|). Such a case was discussed and solved at the end of the 
nineteenth century by Liouville, Appell and others and is presented in textbooks such as 10, [TT| . 
More generally, for polynomial ODEs of the form 

y' = f n (x)y n + h(x)y + f (x) (4) 

Chini |(| |l0) presented at the beginning of the twentieth century a method similar to this mapping into 
separable but through transformations (^|) with q = 0. Chini's method is equivalent to solving (^) by 
determining, when they exist, symmetries (Q) with Q = 0. 

In connection with the above, this work presents a generalization of these methods as an algorithm for 
determining whether or not a first order ODE of arbitrary form|^| 

y' = *(a;, y), 

lr The Abel ODE members of a given class can be mapped between themselves through (0). There are infinitely many 
non-intersecting such Abel ODE classes but only some of them (still infinitely many) admit symmetries of the form (|f]). 
2 R,iccati ODEs are partially excluded from the discussion - see sec.H. 
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belongs to this class admitting linear symmetries (|T|), and, if so, for explicitly finding the symmetry, without 
restrictions to the form of {F, P, Q}. Both the determination of the existence of a solution as well as of the 
symmetry itself are performed without solving any auxiliary differential equations. 

The exposition is organized as follows. In sec.|^, the connection between the symmetries of the form ([!]) 
and linear transformations of the form (||) is analyzed and a solving algorithm for the related ODE class 
is presented. So me e xamples illustrating the type of ODE problem which can be solved using this method 
are shown in sec. 2.2. In secTj, a discussion of Riccati ODEs in terms of their symmetries is made and a 



variant of the method of sec.|2j, to solve a subset of the Riccati problem, is shown. In sec.|J, a discussion and 
some statistics related to the classification of Kamke's first order ODE examples is presented. Finally, the 
conclusions contain general remarks about this work. 



2 Linear transformations and symmetries 

To determine whether or not a given first order ODE has symmetries of the form (Q), following we take 
advantage of the fact that for such a symmetry, the related invariant ODE family can be computed in closed 
form. With the invariant family in hands we then show how one can algorithmically compute infinitesimals 
(Q), when they exist, by quadratures; i.e., using only elementary algebraic manipulation, differentiation, 
integration and constructing inverses of maps. For that purpose, with no loss of generality, we first rewrite 
(0) in terms of {f (x) , p(x) , q(x)} entering (g) as[] 

1 p'y+g' 

7?= — fT (5) 

A direct computation of the finite transformations generated by (|5|) shows that they are linear, of the 
form 

, v P(x)y + q(x) -g(X) 

t = x ' u = Rx) ( ) 

where X is a solution of /(X) — f(x) — a = and a is the (Lie) group parameter. 

Regarding the invariant ODE family related to (|5|), this family can be obtained, for instance, by computing 
differential invariants {lo, I\} of order and 1 related to (g): 

/' 

h=PV + q, h = — : — 7 (7) 

p y + py + q 

and hence, the invariant ODE, given by A(io, I\) — 0, where A is arbitrary, can be conveniently written as 
IiG(Iq) — 1, with arbitrary G, resulting in 

tf=L G {py + q)-t-£y (8) 
P P P 

Due to this connection between linear transformations and symmetries (|l|), the same invariant family (^|) 
can be obtained directly from an autonomous ODE, 

v! = G(u) (9) 

by just changing variables in it using (|^). The solution to (^) can be obtained in the same way, by changing 
variables in the solution to (||): 

rpy+Q i 

f-J W) dz + (h=0 (10) 

In fact (]|) too can be obtained from the symmetry [£ = 1, r\ = 0] of (p|) by changing variables using (||). We 
recall that the knowledge of {£, rj} entering (||) suffices to express the solution of any member of the class 
(||) by quadratures and so it is equivalent to the determination of {/, p, q} entering (10). 

3 In this section we assume /' ^ 0, since, otherwise, (^) would be a first order linear ODE. The limiting case where 
/' — + oo, £ = 0, however, is not excluded. 
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Theorem 1 Consider an ODE y' = ty(x,y) with ^ yyy ^ 0. Both determining whether or not this ODE be- 
longs to the class and, in the positive case, computing the infinitesimals ^) themselves, can be performed 
algorithmically from ^ by quadratures. 

Proof. We start by noting an intrinsic feature of ODEs that are members of the class (||) : the determination 
of p(x) up to a constant factor - say n - suffices to map any member of the class into another one having a 
symmetry of the form 

£ = F(x), V = Q(x) (11) 

This is easily verified by changing variables 

y = — (12) 

up 

in (^), arriving at [£ = l/f, rj = —q'n/f], which is of the form (pi] ). In turn, symmetries of the form (|Tl|), 
when they exist, can be systematically determined as shown in M. In what follows we develop the proof by 
first showing how to map any ODE member of (^|) into one having a symmetry of the form (|TT|), and then, 
for completeness, briefly reviewing how such a symmetry is determined when it exists. 
Now, in view of (|[), for all ODEs y' = ^ members of that class we have that 

*„ = f'G'(py + q)-£ 

* vv = fpG"(py + q), (13) 
*«/» = fp 2 G"'( P y + q). 

Thus let 

A = p^ = -K(py + g) where K=%. (14) 

Vyyy Lf 

Three cases now arise, respectively related to whether A y = 0, A yy = or A yy ^ 0. 
Case Ay = 

In this case, K' = 0, so that K — n for some non-zero constant k, and hence A — n/p. So, from (|l2|), when 
y' = $ is indeed a member of this class, by changing variables using 

y = Au (15) 
the resulting ODE family will have a symmetry of the form (pd|). 

Case A yy = 0, A y ^ 

In this case, K" — so that from (|l4| ) 

A = n 1 (py + q) + k q 
P 

for some constant ko and some non-zero constant K\. Here the necessary condition for y' = to be a 
member of (||) is that the ratio above be linear in y. In such a case, when the ODE is indeed a member of 
this class, by introducing 



u = \n{A) (17) 



the resulting ODE in u will have a symmetry of the form (|1 1|) ; this can be verified straightforwardly by 
performing the change of variables directly in (pf) . 
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Case A yy ^ 
In this case let 

rlii A- n' 

(18) 



_ A yx _ p'y + q' 



Ayy P 

The necessary condition for y' = ^ to be a member of (@) is then that / be linear in y, from where 



p(x) — exp Iydx^j 



(19) 



So, when the ODE belongs to this class, from (|l2|), changing variables y = u/p will lead to an ODE having 
a symmetry of the form ([ll]) . 

Once we have shown how a member of (0) can be mapped into another one having a symmetry of the 
form (|ll|), what remains to be done in the proof of Theorem |l| is to review how that symmetry can be 
obtained by quadratures. 

2.1 Symmetries of the form [£ = F(x), rj = Q{x)\ 

By computing differential invariants, as we did to arrive at (||), the invariant ODE family associated to 
[£ = F{x), r\ = Q(x)] can be written asQ 

y> = $(x, y) = -Jj-? (q(x) +G(y- f ^rdx] ) (20) 



F(x) V V J F ( x ) 

where F, Q and G are arbitrary functions of their arguments. So far we have shown that if an ODE belongs 
to (|8|), then after changing variables as shown for the cases A y = 0, A yy — 0, A yy ^= 0, the resulting ODE 
will be a member of this family (20). 

Now, to determine F and Q, following 0, we first build an expression depending on x and y only through 

G 

K =®l_ = <b_ (21) 
©tt Stt 

where we assume^] <S> yy ^ 0. As explained in 0], the problem then splits into two cases. 
Case /C| ^ / 

In this case, we can obtain the ratio Q{x) / F(x) - only depending on x - by taking 

/C-f t (x) 

The knowledge of this ratio in turn permits the elimination of Q from the determining PDE for the infinites- 
imals, leading to| 

F(x) = Ci e x 7 (23) 

which together with (p2| ) gives the solution we are looking for. The necessary and sufficient conditions for 
the existence of such a symmetry are: 

±(K±) =0 ± ( r*v-r*-** \ (24) 

SyUJ ' dy{ $ + T J" (24) 

4 If Q = or F = then the invariant ODE is separable or linear, so in ( pc| ) and henceforth we assume F ^ 0, Q ^ 0. 
5 If & yy = 0, then ( p^ ) is already a first order linear ODE solvable in terms of quadratures. 
6 For more details see [W 
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Case K\ = I 

Since /Cj = — IC^J-fQ, then when JCj = /, /C§ = / too, so /C = k for some non-zero constant k. Hence, the 
right-hand-side of (^) satisfies: 

^yy 

and so ( p0| ) - the invariant ODE family - is of the form 

y' = $ = A{x) + B(x) e^^ (26) 

where A and -B are arbitrary functions. For a given ODE of this type, A and B can be determined by 
inspection, and the determining PDE for the infinitesimals can be solved directly in terms of A and B as 

— dx 

K 

F(x) = , Q(x) = AF(x) (27) 

2.2 Examples 

1. Consider the first order ODE example number 128 from Kamke's book: 

xy' + ay- f(x)g(x a y) = (28) 

For this ODE, Kamke shows a change of variables mapping the ODE into separable, derived for this particular 
ODE family in 0. Using the algorithm presented in this paper, we tackle this ODE by computing A in ( |l4| ) 
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x a g"' 



so we are in "Case A yy ^ 0". We then proceed with computing / (see (|18j)) arriving at 

j _ ay 

x 

The existence condition that I be linear in y is satisfied, hence, according to (jl9|), 



(29) 



(30) 



p(x) = x a (31) 



Now, changing variables y = — as indicated in (|l2]), (pj|) becomes 



u'=g{u) f -^ (32) 



which is already separable (and thus naturally has a symmetry of the form (|ll|)). The solution to ( |28| ) is 
then obtained by changing variables back in the solution to ([32]), leading to: 

rx a y i 

x a - 1 f(x)dx-j -^-jdz -G\=Q (33) 



2. Let's now discuss an example for which A yy = 



y' = ( x 3 y A + 4 x 4 y 3 + 6 x 5 y 2 +4x 6 y + x 7 ) (x n + 1) - - - 2 (34) 

x 



where n is an arbitrary constant. For this ODE, from ([Ljj), 



A = »±£ (35) 
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so the change of variables here is u — \tl{A) (see (|l7|)), mapping (34) into 



u , = t 7+n ^n + 1 ) e 6« + 



For this ODE, a symmetry of the form ( |11| ) is computed algorithmically (see sec. 2.1) 

1 1 



8 (t n + 1) ' 8t(i n + l) 
Changing variables back directly in the above we arrive at a symmetry for ( |3~i| ) 

1 y + 2x 



S 8(a;" + l) : 
from where an implicit solution to ( |34| ) follows as 

1 



V = 



\x(x n + 1) 



3 x 3 (y + xf 1 + n 
3. As an example of the case in which A y — 0, consider 

(x 2 -l)y 1 
v' = be axy x a + ± '— - - 



= C X 



2 + hi(a;) + c 



(36) 



(37) 



(38) 



(39) 



(40) 



x x 

where a, b and c are arbitrary constants. From (|IJ), A = l/(ax), so that by changing variables as indicated 
m (|l|), © becomes 

(41) 



/ = ut + abt {a+1) e u - -+at{\n(t) + c) 



and this ODE has a symmetry of the form (|1_ 

, 1 

from where, by changing variables back, (^o|) admits the symmetry 



a; ?/ + 1 



(42) 



(43) 



which suffices to integrate (Q) by either using canonical coordinates or computing an integrating factor. 

4. The algorithm presented is applicable to higher degree ODEs too (see Table 1. in sec.||), provided 
that it is possible to solve the given ODE for y' . Consider for instance Kamke's example 394 



( y >f + 2fyy' + gy 2 + (f 2 -g)e- 2 f: f{z)dz 



o 



(44) 



where / = f(x) and g = g(x) are arbitrary functions. For this problem Kamke presents a particular change 
of variables derived in |l2[. To tackle this example using our algorithm we first solve the ODE for y' 



v' = -fy ± 



y 2 - e 



(45) 



Now, by taking any of the two branches of ((45[), for instance the "+" one, and computing the second 
derivative of ( 14 ) , we find 

Ayy = — z ~~2 (46) 



3y 3 e 
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so that I in (|T^) is given by 

i = fy (47) 

from where we compute p using (|l9|), finally arriving at the symmetry 

£ = ^=> v = -^= (48) 

actually admitted by both branches of (fi"5|). 

It is worth mentioning that if on one hand these four examples are straightforward problems for the 
single solving algorithm presented, on the other hand we are not aware of any other algorithm for tackling 
examples like 2 or 3; also, for examples 1 and 4 the changes of variables presented in Kamke are non-obvious 
and presented in connection with different problems [Q, [l2| . Despite the presence of arbitrary functions and 
parameters, both Kamke's examples 128 and 394 are actually particular cases of the class represented by 

(§)■ 

3 Riccati equations 

The case of Riccati type ODEs 

V' = hv 2 + .hv + k (49) 

where fi = fi(x), f% ^ and fo ^ 0, deserves a separate discussion. All Riccati ODEs admit symmetries 
of the form ([!]) and so all of them can be mapped into separable using transformations of the form (||). 
However, it is easy to verify that to find such a transformation requires solving the Riccati ODE itself. The 
algorithm of the previous section - which does not rely in solving auxiliary differential equations - only works 
when ^yyy / - see Theorem 1. The usual approach for solving ( ^9| ) then consists of converting it to a 
linear second order ODE and using the various methods available for this other problem. 

Nonetheless, there are entire subfamilies of (|l9|) for which symmetries of the form ([]]) can be found 
following an approach such as the one presented in the previous section, without using techniques for linear 
second order ODEs. Such an approach is interesting since it enriches the algorithms available for tackling 
the problem and could be of use for solving some linear ODEs by mapping them into Riccati ones as well. 
For the purpose of discussing these cases and without loss of generality we rewrite (||) - the general form of 
the symmetries admitted by Riccati ODEs - redefining /' — > / jp 

p p'y+q' f . n , 
f= ? 1= — (50) 

If now, in (|J) , we redefine /' in the same way and take G as a square mapping depending on two arbitrary 
constants a and 6, 

G = u i— » u 2 + au + b (51) 

we arrive at the form of an arbitrary Riccati ODE - as general as ( fl9l ) - but expressed in terms of these two 
constants a and b and the functions {/, p, q} appearing in its symmetry generator (^p|): 

y ,= f f + (" + >*)f-P' y+ «* + «)«+J>)f-*P (52) 
p p z 

Case p' = 

A first solvable case happens when p' = 0, so that in ([so|) both infinitesimals depend only on x and hence 
the symmetry can be systematically determined - also for Riccati ODEs - as explained in sec. 

Case q' = 



2.1 



A second solvable case happens when, in (|52j), q' = 0, so that the infinitesimals (50) are of the form 



Z = F{x), r, = P(x) y (53) 
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An algorithm for solving such an ODE was presented by Chini ||. In the case of Riccati ODEs, Chini's 
algorithm can be summarized as "to check for the constant character" of the expression]] 

.fof'2 — /0/2 - 2 /0/1/2 , 

X = ± , '— (54) 

(io/2) 3 

where fi are the coefficients of y in (|49|). Whenever X is constant, the problem is systematically solvable in 
terms of quadratures (see for instance Q - p. 303). Concerning ( p3^ ) at q' = 0, a direct computation of X 
confirms that in such a case X is constant. Conversely, another direct computation shows that whenever T 
is constant, the ODE will have a symmetry of the form (p3). To check that, it suffices to solve (p4) for f\ 



and substitute the result into (49); the resulting Riccati ODE will admit the symmetry 




(55) 



which is of the form ( p3| ) . 

We can also see what is the ODE class solved by Chini's algorithm, as well as explain the previous results, 
by noticing that (E3) is an absolute invariant for (B^) under transformations of the form 



that is, of the form 



where 



t = f(x), u=p(x)y 
with g = 0. In fact, ( pi| ) can be written as 

2 



I 



S3 



S2 



(56) 

(57) 
(58) 



/0/2, S3 — /o /2 — /2 /o — 2 /0/1/2, 

are relative invariants of weight 2 and 3 with respect to transformations ( |56| ) Q . 

In summary: Chini's algorithm solves the ODE class generated by changing variables ( |5^ ) in the general 
Riccati ODE ([52]) at q' = 0, all of whose members have X = constant. 

Three additional solvable Riccati families, where the invariant X is non-constant, are obtained by equating 
in d50| ) any two of the three arbitrary functions {/, p, q} with p' ^ and q' ^= 0. 

Case f = p 



When a given Riccati ODE belongs to this family, then by changing variablesj 

u 



y = 



f 



(59) 



in the given ODE and in the general form of its symmetry (|5C|), we see that the resulting ODE in u will 
admit the symmetry 



f = i, v = -q' 



(60) 



which can be determined as explained in sec. 2.1. 
Case q = p 

When a given Riccati ODE belongs to this family, then by changing variables 

y = u-i (61) 

7 This connection between the constant ch ara cter of an expression like built with the coefficients of a polynomial ODE 

of the form (kh, and symmetries of the form (p3h, is valid not only for Riccati ODEs but for ODEs of the form (W) in general. 
8 / is the coefficient of y 2 in the given ODE. 
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in the given ODE and in (^) we see that the resulting ODE in u will admit a symmetry of the form 

£ = j; V = -j u (62) 

that is, infinitesimals of the form (|5^), and hence the ODE will be solvable using Chini's method. 
Case / = q 



From (52), the Riccati family corresponding to this case is given by 



y , = tf + (f(" + *f)-P , ) y+ «° + f)f+J>)f-f'P ( 63 ) 
p p z 

From (^0|), this ODE family admits the symmetry 

e p P'y+f fR A\ 

7' 11 = 7 — 

We haven't found an obvious transformation of the form y = P u + Q to map this symmetry to one of the 
forms (|ll|) or (|5^) . A possible approach would then be to directly set up the determining PDE 

*k + (vy - £') {hv 2 + hv + Jo) - £ (/2'y 2 + h'y + fo) - v (2 hv + h) = 

for the coefficients ^ and 77 of the infinitesimal symmetry generator £^ + 77^- of an arbitrary Riccati ODE 



(49). Then take £ and rj as given by Q64| ) and run a differential elimination process solving for p. That 
approach works but the resulting symbolic expressions are large enough to become untractable even with 
simple examples. 

An alternative approach leading to more tractable expressions is based on using the information we have 



in (53) concerning the existence of two constants a and b. So, by equating the coefficients of (63) with those 
of (49) we arrive at the system 

fa ~ f = 0, (65) 

, / 2 (a + 2/ 2 )-p' n . / 2 (/ 2 « + /2 2 + b)-/ 2 p n 

h = u, Jo 5 — u 

p p z 

This system can be solved for a and b, so that when a Riccati equation belongs to this family the following 
two expressions formed with its coefficients fi will be constants: 



hp-2f 2 2 +p' 

f 1 n — n ) -\- fn n 

(66) 



h 

hp 2 +f2(f2 2 -hp- P ')+hp 



h 

At this point, however, we cannot verify the constant character of these expressions because p is still unknown. 
An expression for p can be obtained by computing the integrability conditions for /a and p implied by (|6q) : 

, f/ 2 ') 2 +p/o/ 2 ' 

h - hh -2/ p'-p/o + V ' ; , (67) 



2 



2hh-p' h-ph + hP ' + t Phh (68) 

h 
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Using (|67]) to eliminate p' from (|6^) we arrive at a solution for p: 

V = (69) 



h I 3s 2 / 2 s 2 - 2s 2 2 / 2 + ^(s 2 - s 4 ) s 2 - 8s 2 d - 2 (a a J + 2s 3 2 j h 

si (2s 2 s 3 ' - 3s 3 s 2 ') 

where 

2s 2 s 3 ' - 3s 3 s 2 ' + 3s 3 2 

s 4 = = 70 

2 s 2 

is the relative invariant of weight 4 for Riccati ODEs with respect to transformations of the form (p6|) H . 



We note that, from (p7|), when the denominator of (|69j) is zero the ODE has a constant invariant X (54) and 
so it is already solvable using Chini's algorithm. 

In summary, a strategy - not relying on solving second order linear ODEs - for finding linear symmetries 
of the form (Q) for Riccati equations could consist of 



1. Check if the ODE has a symmetry of the form [£ = F(x), r\ = Q(x)] (algorithm of H - see sec. 2.1); or 
[£ = J-{x), Tj = V(x) u] (Chini's algorithm); 

2. Check if the ODE belongs to one of the two families "/ = p" or "p = q" by using the transformations 
j59| ) and (^lj) and re-entering the previous step; 

3. Check if the ODE belongs to the family "/ = q" (|63|); for that purpose: 

(a) compute the invariants (p8[); 

(b) use these invariants to compute p using (^9|) ; 

(c) plug the resulting p into ( p6| ) and verify if the two right-hand-sides are constant - if so, the ODE 
admits the symmetry ( p4| ) 

It is our belief that with the development of computer algebra software and faster computers, these type 
of algorithms for Riccati subclasses (there are of course many other possibilities) will become each day more 
relevant, also as an alternative for tackling second order linear ODEs. 



4 Classification of Kamke's examples and discussion 

We have prepared a computer algebra prototype of the algorithms presented in sec.|^ and sec.|| using the 
Maple system. We then used this prototype to analyze the set of Kamke's 576 first order ODE examples. 
This type of computational activity interesting. In the first place, Kamke's book contains a fair collection of 
examples arising in applications. In the second place, this book collects many, perhaps most, of the solving 
algorithms available in the literature for first order ODEs. So, generally speaking, such an analysis conveys 
a reasonable evaluation of the range of application and novelty of a new ODE solving algorithm. 

In order to perform the classification, we first excluded from the 576 Kamkc examples all those for which 
a solution is not shown and we were not able to determine it by other means^j - most of them just because 
the ODE is too general. So our testing arena starts with 552 ODEs. 

Then, we know that all ODEs of type separable, linear, homogeneous, Bernoulli, Riccati and Abel with 
constant invariant)^], that is, 372 of Kamke's examples, have symmetries of the form ([!]) and then belong to 
the ODE class discussed in this paper. So the first thing we wanted to know is how many of the remaining 
552 - 372 = 180 examples admit linear symmetries of the form ([!])? Is there any other classification known 
for these 180 examples? The information for answering these and related questions is summarized in this 
table: 

9 The numbers of the Kamke examples we excluded in this way are: 47, 48, 50, 55, 56, 74, 79, 82, 202, 205, 206, 219, 234, 
235, 237, 265, 250, 253, 269, 331, 370, 461, 503 and 576. 

10 For an enumeration of Kamke's examples of Abel type with constant invariant see [fej and concerning other classes sec a\. 
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Class 


First degree in y'\ 8> 


S Higher degree in y': 92 


Total: 180 ODEs 


K = F, V = Py + Q] 


20 


37 


57 


Abel (non-constant invariant) 


15 





15 


Clairaut 





15 


15 


d'Alembert 


2 


21 


23 


Unknown 


33 


21 


54 



Table 1. Classification of 180 non "linear, separable, Bernoulli, Riccati or Abel c.i." Kamke's examples. 



From these numbers, some first conclusions can be drawn. First, 372 + 57 = 429 ODEs, that is, 78 % of 
Kamke's 552 solvable examples, have linear symmetries of the form (Q). Also Table 1. shows that in Kamke 
- even among these particular 180 ODEs which exclude the easy ones - there are more examples having 
symmetries of the form (|l|) which can be systematically determined as explained in sec.|^, than examples of 
Abel (non-constant invariant), Clairaut and d'Alembert types all together^. 

In the second place, if we discard the 61 examples of Riccati type found in Kamke, the solvable set is 
reduced to 552 - 61 = 491 ODEs, and from this set, 429 - 61 = 368 ODEs, that is, 75 % of these 491 can be 
solved algorithmically as shown in sec.|^. 

Moreover, a classification of Kamke's examples of Riccati type according to sec.|| shows: 



Class [£_=F\ V = Q] [g = F, r] = Py] "Two of {f,p, q} are equal" Total of ODEs 

Riccati 7/61 22/61 2/61 31/61 

Table 2. Classification according to sec.^j of the 61 Riccati ODE examples of Kamke's book. 



So one half of these Riccati examples is still solvable using the algorithms described in sec.g, and so without 
mapping the problem into a second order linear ODE nor having to solve auxiliary differential equations. 



5 Conclusions 

In this work we presented an algorithm for solving first order ODEs, consisting of determining symmetries of 
the form (jl]). From the discussions of sec.|[ to these symmetries one can associate an ODE class - represented 
by (||) - which embraces all first order ODEs mappable into separable ones through linear transformations. 
From the numbers of sec. [|, this class appears to us as the widest first order ODE class we are aware of, all 
of whose members are algorithmically solvable as shown in sec. || or sec. [|, or mappable into second order 
linear ODEs when the ODE is of Riccati type and the methods in sec. | don't cover the caseQ. 

The algorithms presented also do not require solving additional differential equations nor do they rely 
on the ODE member of the class being algebraic (i.e.: rational in y and its derivatives) or on restrictions to 
the function fields - the only requirement on the functions entering the infinitesimals (|l|) are those implied 
by the fact that these infinitesimals do generate a Lie group of transformations. 

Concerning other related works we are aware of, the method presented in Jl3| for solving Abel ODEs 
with constant invariant is a particular case of the one presented here in that those ODEs are the subclass of 
(||) of type Abel. In the same line, the method by Chini [|| - a generalization of the method for Abel ODEs 
with constant invariant which solves more general ODEs of the form (^) - is also a particular case in that 
([I|) is a very restricted subclass of (^). Also, the class solvable through Chini's method is equivalent to a 
separable ODE only through transformations (Q) with q = 0. In this sense, the algorithm presented in sec. ^| 
generalizes both the one discussed in |lj| and the one presented in Q . 

The fact that this class (||) is algorithmically solvable turns this classification relevant for modern com- 
puter algebra implementations too. As shown in sec.|4|, taking as framework for instance Kamke's examples, 

11 We note there is no intersection between these classifications: all of Abel (non-constant invariant) Clairaut and d'Alembert 
equations do not have symmetries of the form (hi). 

12 In fact the solving of the equivalence problem for Riccati ODEs with respect to linear transformations is equivalent to 
solving the whole Riccati class; in turn a too general problem. 
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78% belong to this ODE-class. Even after discarding Riccati ODEs, 75% of the remaining Kamke examples 
belong to this class (||) and so are solvable by the algorithm presented in sec.|[ This algorithm actually 
solves many ODE families not solved in the presently available computer algebra systems (CAS). For in- 



stance, from the 4 examples shown in sec. 2.2, three of them cannot be solved by Maple 6 or Mathematica 4 
- the last versions of these major CAS. Even concerning the Riccati families presented in sec.^J, two of them 
(cases "/ = g" and "p = q") are also not solved by these two CAS, which base their strategy in mapping 
Riccati ODEs into linear second order ones. 
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